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COLLISIONLESS  TEARING  INSTABILITY  OF  A  BI-MAXWELLIAN 
NEUTRAL  SHEET  WITH  EXACT  PARTICLE  ORBITS 


I  Introduction 

The  collisionless  tearing  mode  in  a  neutral  sheet^“®  has  received  considerable  at¬ 
tention  in  the  past  two  decades®" A  neutral  sheet  can  be  described  by  a  one¬ 
dimensional  magnetic  field  B(r)  =  Bx(:)ic  such  that  Bx  is  odd  in  r,  the  coordinate 
normal  to  the  plane  of  the  sheet.  A  well-known  model  is  the  Harris  eciuilibrium^^ 
which  is  established  by  balancing  J  x  B  and  Vp,  where  J(2)  =  Jy[z)y  is  in  the  plane 
of  the  neutral  sheet.  In  its  classic  form^®  the  instability  is  due  to  the  finite  mass 
of  the  current  carriers,  allowing  the  current  sheet  to  filament  (k  —  kx'k)  and  attain 
a  lower  magnetic  energy  configuration.  For  this  process,  the  existence  of  regions 
about  c  =  0  (the  null  plane)  in  which  electrons  and  ions  are  essentially  unmagne¬ 
tized  is  critical.  As  the  magnetic  field  is  perturbed,  an  electric  field  is  induced  in 
the  y  direction,  accelerating  the  unmagnetized  electrons  and  ions.  In  our  model,  we 
will  ignore  the  scalar  potential.  This  instability  is  potentially  important  for  mag¬ 
netic  reconnection  processes  that  may  occur  in  a  number  of  physical  systems  such 
as  the  magnetotail*^'*®  and  in  the  day-side  magnetopause^^  *®.  In  realistic  physical 
systenis,  the  magnetic  field  is  generally  not  one  dimensional.  For  example,  the  mag¬ 
netotail  would  have  a  small  B-  component  and  the  magnetopause  would  have  a  By 
component  comparable  to  the  asymptotic  Bx  component.  However,  the  instability 
occurring  in  a  neutral  sheet  may  be  prototypical,  elucidating  the  essential  physics 
common  to  various  configurations. 

The  majority  of  the  previous  papers  on  the  collisionless  tearing  mode  have  used 
Maxwellian  equilibrium  distribution  functions,  for  which  the  mode  grows  only  slowly 
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(-y/wc,  "Cl).  Applied  to  the  earth’s  niagnetotail,  for  example,  this  gives  e-folding 
times  which  are  significantly  longer  than  the  relevant  time  scales  for  possible  I’econ- 
nection  except  in  extreme  situations.  For  example,  using  the  Harris  equilibrium, 
Coppi  et  al.®,  estimated  7  =sl/.5  sec“*  -  1/1.5  sec“'  with  pi/ 8  %  1,  where  pi  is  the 
asymptotic  gyroradius,  and  8  is  the  magnetic  scale  length,  which  requires  that  the 
global  bulk  plasma  flow  across  the  tail  be  comparable  to  the  thermal  velocities.  Such 
high-speed  bulk  flows  are  not  necessarily  observed  prior  to  reconnection.  However, 
it  has  been  know’n  that  an  electron  temperature  anisotropy  can  significantly  affect 
the  tearing  instability  growth  rate.  Forslund’®  showed  that  the  mode  was  stabi¬ 
lized  in  equilibria  with  (1  -  Tei/Te||)  >  P^I8.  However,  Laval  and  Pellat*'  later 
showed  that  these  equilibria  were  unstable  to  a  quickly  growing  obliquely  propa¬ 
gating  mode.  A  weak  anisotropy  >  7’e||  greatly  increases  the  growth  rate^®.  In 
addition,  an  anisotropic  neutral  sheet  should  be  unstable  to  the  Weible  instability*’’. 
Thus,  any  electron  anisotropy  in  the  magnetotail  plasma  is  likely  to  be  non-linearly 
isotropized  by  these  and  other  instabilities  on  much  faster  time  scales  than  the  typ¬ 
ical  time  scale  associated  with  reconnection  in  the  magnetosphere.  More  recently, 
Chen  and  Palmadesso*  investigated  the  dependence  of  the  growth  rate  on  the  ion 
temperature  anisotnipy.  They  found  that  an  ion  temperature  anisotropy  could  also 
enhance  the  tearing  mode  growth  rate.  For  T,i_/Ti\\  ~  1.2  —  1..5,  the  growth 

rate  was  fottnd  to  be  one  to  two  orders  of  magnitude  greater  than  in  the  isotropic 
case.  This  was  later  confirmed  by  simulation  results***.  For  small  to  moderate  ion 
anisotropy  r',i/T’,||  1..5,  the  e-folding  time  can  be  reduced  to  a  few  minutes  for  the 

earth’s  magnetotail  parameters,  which  is  favorable  in  comparison  with  the  observed 


2 


substorm  onset  time  scales  of  tens  of  minutes.  This  suggests  that  non-Maxwellian 
features  in  the  ion  distribution  function  rather  than  in  the  electron  distribution 
are  of  interest  for  magnetospheric  processes.  We  do  not  expect  a  significant  ion 
temperature  anisotropy  to  be  present  in  the  quiescent  magnetotaiP®'^^’^^.  However, 
in  a  magnetotail-like  geometry  phase  space  is  partitioned  into  disjoint  subregions 
containing  transient,  chaotic,  and  integrable  orbits^^.  Information  of  disturbances 
far  from  the  current  sheet  would  reach  these  subregions  in  different  characteristic 
times,  leading  perhaps  to  non-Maxwellian  features  in  the  distribution  function,  a 
process  referred  to  as  “differential  memory”^^  These  non-Maxwellian  features  can 
persist  on  time  scales  of  orbit  diffusion  across  phase  regions,  a  longer  time  than  the 
e-folding  time  of  the  ion  anisotropic  tearing  instability. 

In  works  previous  to  and  including  ref.  [l]  the  complicated  particle  orbits  of 
the  neutral  sheet  geometry  have  only  been  accounted  for  in  approximate  ways. 
One  method  of  approximating  the  orbits  is  a  global  energy  principle®  ' computed 
to  first  order  in  the  ratio  of  the  Larmor  radius  to  the  magnetic  scale  length  (the 
electron  Larmor  radius  for  an  isotropic  plasma).  This  method  has  the  advantage 
that  a  particular  form  for  the  functional  dependence  of  the  perturbation  on  z  is 
not  assumed  but  the  disadvantage  that  only  a  stability  criterion  may  be  found.  A 
secotid  method  uses  both  approximate  orbits  and  assumes  a  particular  form  for  the 
eigenfunction.  The  most  common  example  of  the  latter  method  is  the  two  region 
approximation*^'^  ’”'^®.  In  this  approximation  all  orbits  within  a  nearby  region  of 
the  field  null  are  taken  to  be  straight  line  orbits,  and  the  eigenfunction  of  the 
field  perturbation  is  assumed  to  be  constant  within  that  region  (the  “constant  psi” 
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approximation).  Outside  of  this  region,  the  orbits  are  taken  to  be  the  usual  drifting 
cyclotron  orbits.  Provided  the  ion  distribution  is  isotropic  and  (pe/Pt)  “C  1,  where 
Pe  and  Pi  are  the  Larmor  radii  of  the  electrons  and  ions,  the  ion  contribution  to  the 
non-adiabatic  perturbed  current  is  small®.  Thus,  attention  need  only  be  paid  to  the 
crossing  electron  orbits. 

An  interesting  facet  of  the  work  in  ref.  [1]  was  the  identification  of  an  ion- 
intermediate  region  within  which  the  ion  orbits  were  approximated  by  straight  lines 
but  the  electron  orbits  were  approximated  by  drifting  guiding  center  orbits  (the 
“3-region”  approximation).  This  ion-intermediate  region  proved  to  be  important  in 
the  ion  anisotropic  case;  if  the  region  was  included  the  tearing  mode  growth  rate 
increased  by  nearly  an  order  of  magnitude  over  a  similar  calculation  where  the  ion 
contribution  to  the  current  outside  of  the  electron  inner  region  was  taken  to  be 
negligible  (the  usual  “2-region”  approximation).  It  was  therefore  expected  that  the 
straight-line  approximation  for  orbits  within  this  large  ion-intermediate  region,  as 
well  as  the  “constant  psi”  approximation  could  be  limiting  for  all  but  the  smallest 
values  of  p,  in  the  presence  of  an  ion  anisotropy.  Recently,  Chen  and  Lee*®  have 
developed  a  method  in  which  growth  rates  could  be  calculated  taking  into  account 
all  of  the  orbits  exactly.  This  method  removed  the  limitation  in  applicability  arising 
from  the  constant-^’  approximation  and  from  approximate  orbits.  They  applied  the 
treatment  to  a  delta  function  particle  distribution  in  perpendicular  energy.  This 
choice  gave  rise  to  a  highly  non- Maxwellian  distribution  function  with  a  very  fast 
growth  rate,  i.e.  approaching  the  ion-cyclotron  frequency.  The  eigenfunction  was 
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highly  structured  within  the  electron  inner  region,  so  that  the  constant-^’  approxi¬ 
mation  would  have  been  invalid.  However,  such  a  specialized  distribution  function  is 
not  likely  to  occur  in  physical  systems,  such  as  the  magnetosphere;  quantitative  re¬ 
sults  recjuire  a  more  broadly  applicable  distribution  function  such  as  a  bi-Maxwellian 
distribution  function.  In  addition,  their  choice  has  not  been  previously  analyzed  by 
approximate  linear  theories,  or  by  simulation,  and  there  is  no  isotropic  limit. 

In  this  paper,  we  apply  the  integro-differential  method  of  Ref.  [26]  to  the  bi- 
Maxwellian  system.  In  addition  to  being  more  physical  for  the  magnetotail,  this 
model  has  the  advantage  that  the  equilibrium  quantity  is  continuously  vari¬ 

able  so  that  it  would  be  a  better  basis  for  a  more  complete  study.  We  have  solved 
the  full  integro-differential  equation  for  the  bi-Maxwellian  neutral  sheet  with  no 
restrictions  on  T',i/r,||  or  Pi/S.  The  solution  shows  that  for  Til/7’,11  >  1  the  growth 
is  significantly  enhanced  and  the  maximum  growth  wavenumber  k  is  increased  in 
comparison  with  the  isotropic  case,  in  qualitative  agreement  with  an  earlier  study 
(ref.  [1]  which  used  the  3-region  approximation).  However,  the  exact  result  re¬ 
veals  that  the  3-region  approximation  underestimates  the  growth  rate  by  a  factor 
of  two  to  five.  This  is  found  to  be  caused  by  incorrect  classification  of  crossing  and 
non-crossing  orbits  and  by  the  constant-V’  approximation,  which  is  taken  in  many 
approximate  theories. 

Although  the  integro-differential  formalism  can  provide  the  exact  answer,  it  is 
generally  nontrivial  to  use.  In  particular,  for  complex  systems,  the  orbit  contri¬ 
butions  for  each  matrix  element  may  have  to  be  evaluated  repeatedly  in  order  to 
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determine  the  dispersion  relations,  which  may  make  the  actual  computation  pro¬ 
hibitive.  In  this  paper,  we  will  describe  a  new  approximate  technique  which  is  a 
simplified  application  of  the  full  integro-differential  treatment  and  which  properly 
accounts  for  the  essential  features  of  orbit  contributions.  In  particular,  contributions 
from  different  types  of  orbits  are  properly  taken  into  account  and  no  assumption 
on  V’  is  necessary.  This  approximation  allows  much  more  efficient  computation  of 
matrix  elements  with  mitiimal  loss  of  overall  accuracy.  In  order  to  fully  understand 
the  validity  and  advantages  of  this  simplified  integro-differential  treatment,  we  will 
describe  the  the  new  approximation  and  its  physical  basis.  In  this  connection,  we 
will  use  the  three-region  approximation  to  illustrate  the  essential  physics. 

This  paper  will  consist  of  two  major  components.  In  the  first  part  (Secs.  II 
and  III),  we  will  present  the  formulation  and  solution  of  the  full  integro-differential 
equation.  In  the  second  part  (Secs.  IV  and  V),  we  will  discuss  the  reasons  why 
previous  approximations  break  down  and  describe  a  simplified  integro-differential 
method  which  does  not  suffer  from  similar  errors. 
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II  Formulation 


A  Bi-Maxwellian,  Neutral  Sheet  Equilibrium 

A  well  known  equilibrium  exhibiting  a  neutral  sheet  is  the  Harris'*  equilibrium.  For 
each  species  we  choose  a  modified  Harris  equilibrium  which  includes  a  temperature 
anisotropy 


(1) 


where 


exp(  — ) 

F=  ,  -  VP.j/rj.  (2) 

nt'T^ 

Here  =  mvy  +  {q/c)Ay(z),  Hi  ^  mv\/2,  H\^  -  mr||/2,  Ti  =  mvl^/2,  7||  = 
Tnv^/2^  and  V'  is  an  adjustable  parameter  representing  the  average  velocity  of 
the  species  under  consideration.  The  subscripts  “1”  and  “||”  refer  to  components 
perpendicular  and  parallel  to  the  equilibrium  magnetic  field 


B  —  Bo  ♦f>.'ih(r//')e. 


(3) 


The  vector  potential  corresponding  to  this  field  is 


Ay{z)  Bob  In 


cosh(c/^ ) 
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are  a  result  of  the  quasi-neutrality  condition  and  Ampere’s  Law.  Here  p  is  the  gy- 
roradius  in  the  asymptotic  field,  p  =  iYi/w^c  where  =  qB^/mc.  As  a  consequence 
of  these  relations,  w’e  find  the  pressure  profile 

B  The  Fundamental  Integro-Differential  Equation 


In  this  section  we  will  briefly  review  the  integro-diflFerential  equation  for  complete¬ 
ness’  sake.  The  derivation  is  essentially  the  same  as  in^®,  except  that  we  have  the 
simplification  that  dfo/dHw  = 

As  in  the  previous  derivation,  we  take  for  the  perturbed  vector  potential  the 
usual  tearing  perturbation 


(8) 


The  perturbed  y-directional  current  is  found  to  be 


c 
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u;  V 

T„l. 

fioV  "j 

and  Ampere's  Law  may  be  used  to  give  a  single  equation  for  V’ 


,  2  / 

k  li’  -  ^  2  -  ^  j,y 


(10) 


where  Sy  contains  the  integration  over  unperturbed  orbits 
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In  this  configuration,  particle  motion  is  regular  with  a  period  that  depends  upon 
the  constants  of  motion,  T{H  Py).  Thus  the  quantity  Vy(z')il'(z')  may  be  expanded 
as  a  Fourier  series  in  t' 


OO 

n  =  0 


inOf  t*  -t) 


where  Q  =  27r/T.  The  ^•n’s  are  the  Fourier  coefficients  at  the  current  time  f,  that 


IS 


n=0 


These  may  be  found  by  an  integration  starting  at  the  current  time  and  extending 
over  one  period 


1 


-t) 


Hence,  we  may  write  the  identity 

We  may  now  complete  the  /'  integration  in  Sy  to  find 

1  rl  +  T 


-t) 


— u; 

—  7?f2  -  kvil  T 


1  /•'+^ 

J, 


(’  (r") 


Thus  we  have  re[)laced  the  integration  over  an  infinite  number  of  orbits  with  a  sum 
over  an  infinite  number  of  harmonics  of  the  orbital  period.  Physically,  we  know  that 
the  higher  harmonics  n  ^  0  correspond  to  resonances  between  the  orbital  motion 
and  the  oscillation  of  the  eigenfunction.  Resonances  with  the  higher  harmonics  may 
be  discarded  if  the  frequenc_,  of  the  fields  is  much  less  than  the  frequency  of  the 
orbit,  u)  fl  for  all  Py  and  //]  .  When  the  orbits  correspond  to  cyclotron  motion. 
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this  is  identical  to  the  low  frequency  approximation.  In  this  formulation,  we  will 
assume  that  all  contributions  from  resonances  with  higher  harmonics  are  negligible, 
and  therefore  we  ignore  all  of  the  terms  but  n  =  0. 


oj  -  fcrii  T  Jt  L  V"  ^  ^ 


In  the  last  equation  we  have  reduced  the  double  prime  to  a  single  prime  for  conve¬ 
nience.  Our  manipulations  upon  Sy  have  separated  the  parallel  direction  from  the 
perpendicular  direction.  Now,  we  may  complete  the  parallel  velocity  integral  in  the 
perturbed  current 

,  Til  +  JeX  [ 

+  [{Z({)+  (l-  p^)(l+?^(0)]'?(-)}.  (13) 


r  n  V  1  /•^  +  T  r  (  ^  ) 

Q(z)  I  ^  I  d/'[  U'{z')\  (14) 

^  =  u;/fcr||,  and  Z  is  the  plasma  dispersion  function.  The  integration  over  one 
orbital  period  may  be  replaced  by  an  integration  along  the  path  of  the  orbit.  Thus 


^  f  J2  '’y  r.  1  /  (^’y(^') 


The  period  is  given  by 
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C  Method  of  Solution 


The  general  numerical  method  used  here  is  a  reduction  of  the  integro-differential 
equation  to  a  matrix  equation,  and  then  a  numerical  solution  of  the  eigenvalue  equa¬ 
tion  for  the  growth  rate.  Such  a  method  was  used  in  ref.  [26]  so  we  will  only  tovich 
on  this  briefly.  For  further  details  see  refs.  [26,27].  In  our  case,  however,  we  found 
that  we  could  reduce  the  computation,  time  for  each  matrix  element  significantly  by 
changing  the  order  of  integration,  and  we  will  discuss  this  change  in  some  detail. 

As  in  ref.  [26],  the  numerical  solution  of  Eq.  (10)  with  Eq.  (13)  is  found  by 
representing  the  eigenfunction  4'(z)  as  a  sum  over  basis  functions 

=  Y^ar,(i>r,{z). 

n 

Multiplying  Eq.  (10)  by  a  basis  function,  (f>m  a»d  integrating  over  all  z  we  obtain 


0  =  ^  dn^r 


(16) 


where  the  matrix  elements 


and 


^  mv  ^  ^  4^  n  4^  m 

^  _ 

n 

h  T±e 

—  J* 

)  1  Sv 

d<f>n  d<i>„ 

dz  dz 


^  j  dzsec\\^(z/S)(f)„<i>r, 


^^(0+(l-^)(l+^^(^))]A.n  (17) 


jr  ) 


(18) 


We  solve  Eq.  (16)  numerically  by  first  solving  the  eigenvalue  equation 


^  ^  d,x  G ,■ 
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and  then  searching  for  a  frequency  root  of  the  dispersion  relation 


A(a;,A*)=0.  (19) 

In  order  to  find  the  parameters  upon  which  A'„„  depends,  we  find  it  convenient 
to  normalize  our  variables.  When  occurring  within  the  integration  of  A'mn?  time 
will  be  scaled  to  the  cyclotron  period  in  the  asymptotic  field.  Thus  for  the  orbital 
period, 

f  =  u;,T. 

Distance  will  be  scaled  throughout  to  the  magnetic  scale  length  S.  Thus, 

z  =  z/S  and  k  =  kS. 

The  velocities  are  then  scaled  by 

V  —  v/ijJcS. 

Note  that  the  quantity  Vy/V  becomes  Vy(S/p)^-  We  also  scale  the  constants  of 
motion  H±  and  Py  so  that 

h  - 

and 

P  = 

are  only  dependent  upon  i'y,  v.,  and  z.  The  quantity  K„n  becomes 

Krr,„  -  -{fi/pfexp{-{p/Sf]K„„ 

TT 
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where 


and 


Kmn  =  j  dz4>m{z)  j  d^VVyF  \  j  |f,p)|  ) 


F  =  exp[-(6/p)^/i  +  2p\. 


Thus  Kmn  is  a  function  of  p/f)  only.  The  matrix  elements  A'm„  must  be  calculated 
for  each  species  using  the  appropriate  gyroradius  p,  or  pe-  However,  they  need 
not  be  calculated  anew  for  a  change  in  any  other  parameter  such  as  frequency, 
wavenumber,  or  temperature  anisotropy.  This  property  may  not  be  preserved  in 
other  geometries. 


We  now  restrict  the  problem  to  a  consideration  of  only  symmetric  solutions. 
Thus,  we  only  need  perform  the  i  and  z'  integrations  over  contributions  in  the  first 
quadrant.  Therefore,  we  write,  after  a  change  in  the  order  of  integration 


where 

1  p  <  —y/h 

2  y/h  >  p  >  —y/h. 

The  turning  points  and  Z2  are  determined  by  =  0  and  are  given  by 


giVyh)  =  { 


Arccoshjexp 


-y/h  -  pj| 


p  <  -y/h 
y/h  >  p>  —y/h 
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and 


-2 


Arccos)i|exp|\/^  HI 


Note  that  the  value  p  =  —  s/h  unambiguously  distinguishes  between  crossing  (  —  \/h  < 
p  <  \/h)  and  non-crossing  (p  <  ~\/h)  orbits. 


The  factor  g{Pih)  is  due  to  the  distinction  between  crossing  and  non-crossing 
orbits.  This  can  be  explained  as  follows.  For  non-crossing  orbits,  Fig.  1  shows  the 
regions  of  integration  in  the  z  —  z'  plane  for  two  cases  of  the  n  omentum  and  a 
sample  value  of  the  energy.  The  gray  shaded  regions  are  the  areas  of  integration  for 
the  value  p  =  —\/h  —  e  as  e  0.  For  e  sufficiently  small,  non-crossing  orbits  can 
come  arbitrarily  close  to  the  null  plane.  The  cross-hatched  regions  are  the  areas  of 
integration  for  a  smaller  value  of  the  momentum  (p  <  —\/h)  such  that  the  orbital 
turning  points  are  at  £  =  Zi  and  Z2.  Further  decreases  in  the  momentum  p  — >  — oo 
lead  to  farther  separation  in  the  two  areas  of  integration,  the  top  area  tending  to 
z  =  z'  =  oo  and  the  bottom  area  to  z  =  z'  =  —oo.  There  is  no  region  of  integration 
in  the  second  and  fourth  quadrants  (zz'  <  0)  because  no  non-crossing  orbit  passes 
through  two  planes  with  z-coordinate  z  and  £'  that  are  on  different  sides  of  the 
null  plane.  For  crossing  orbits.  Fig.  2  shows  regions  of  integration.  The  interior  of 
the  larger  square  (the  gray  shaded  area  plus  the  cross-hatched  area)  represents  the 
area  of  integration  with  the  largest  value  of  momentum  corresponding  to  a  crossing 
orbit,  p  =  —\/h.  Further  increases  in  the  momentum  lead  to  decreased  areas  of 
integration,  such  as  the  the  cross-hatched  area  shown,  until  the  area  shrinks  to  zero 
for  p  =  \/h,  which  is  the  largest  value  of  the  momentum  kinematically  allowed. 
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Imposing  s5Mumetry  in  f,  only  the  quadrant  with  both  z  and  z'  greater  than  zero 
need  be  included  in  the  integration,  and  therefore  we  simply  must  multiply  the 
crossing  orbit  contribution  by  a  factor  of  four  and  the  non-crossing  orbit  contribution 
by  a  factor  of  two.  The  crossing  orbit  contribution  is  thus  weighted  by  an  additional 
factor  of  two.  For  odd  modes  (in  £),  the  factor  g{p,  h)  is  unity  for  non-crossing  orbits 
and  zero  for  crossing  orbits. 

Equation  (20)  has  three  desirable  properties:  it  is  symmetric  under  interchange 
of  m  and  n,  neither  of  the  coordinates  z  and  z'  is  preferred,  and  for  m  =  n,  the 
result,  Kmni  is  positive  definite.  Numerically,  Eq.  (20)  only  requires  three  nested 
integrations  since  the  z  and  z'  integrations  may  be  performed  independently.  This 
represents  a  substantial  savings  in  computer  time  and  accuracy  over  the  form  used 
in  ref.  [26j  which  requires  four  nested  integrations. 

D  Implementation 


The  basis  functions  which  we  have  chosen  are  the  chapeau  functions. 


f  ( “  — l)/('^m  ~m-l) 

(•"rTT  +  l  *')/(''rn  +  l  ~m  ) 

0 


•^m  — I  _ 

^trt  ^  '^m  +  1 

otherwise 


except  for  the  last  basis  function  at  zyv  >  1 


<i>N 


(z  -  zn_i)/{zn  -  z^r-i) 
0 


^N-l  <  ~  S 

otherwise 


and  the  basis  function  at  the  neutral  plane. 


0  S  -  ■Ji 
otherwise. 
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The  primary  advantage  offered  by  these  chapeau  functions  is  that  the  cen¬ 
ter  points  need  not  be  regularly  spaced.  Thus  in  the  asymptotic  region  where 
V’  ~  exp(  — A;£)  we  may  place  the  basis  function  centers  far  apart.  We  tested  the 
sensitivity  of  our  code  to  the  number  of  basis  functions  by  doubling  this  number  to 
the  present  value  (100)  using  a  mode  with  a  slow  exponential  decay  in  the  asymp¬ 
totic  region  {k  —  0.1)  compared  to  the  fastest  growing  mode.  Such  an  increase 
changed  the  computed  frequency  by  less  than  !j%.  This  suggests  that  the  code  is 
capable  of  finding  very  extended  eigenfunctions  with  good  accuracy. 

Actually,  the  worst  accuracy  limitation  arose  in  the  matrix  solver  (single  pre¬ 
cision  on  a  32-bit  computer),  which  requires  some  diagonality  to  solve  very  large 
matrices.  For  large  pi,  however,  the  large  ion  orbits  make  the  matrix  A'^n  very 
non-diagonal,  and  therefore  G'„„  becomes  non-diagonal.  Since,  in  the  routine  we 
used,  the  accuracy  of  the  computed  eigenvalues  is  decreased  as  the  matrix  becomes 
less  diagonal,  the  frequencies  computed  by  our  method  are  less  accurate  for  larger 
values  of  pi/ 6.  In  practice,  this  led  to  an  inaccuracy  of  around  2%  in  computed 
values  of  the  frequency  for  pi/fi  —  1/2. 
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Ill  Results 


We  have  solved  the  full  integro-differential  equation  Eq.  (10)  by  solving  the  equiv¬ 
alent  matrix  equation  Eq.  (16)  using  exact  equilibrium  orbits.  In  Fig.  3,  the  solid 
line  is  the  solution  of  the  dispersion  relation  [Eq.  (19)]  for  the  isotropic  temperature 
case  (r,i/r,||  =  1)  with  pn'S  —  1/10.  The  dashed  line  is  the  3-region  result  for  the 
same  parameters.  In  this  isotropic  case,  the  3-region  approximation  gives  nearly 
the  same  residt  as  the  classical  2-region  approximation  as  shown  in  Ref.  [1],  and  we 
see  that  the  previous  approximation  agrees  well  with  the  exact  solution.  However, 
it  has  been  shown  that  even  small  to  moderate  deviations  from  exact  ion  temper¬ 
ature  isotropy  can  cause  the  2-region  approximation  to  break  down;  inclusion  of 
the  ion-intermediate  region  can  lead  to  a  one  to  two  order  of  magnitude  increase  in 
growth  rate’  In  Fig.  4,  we  show  the  case  for  Tix/Tlii  =  1.2  with  p,,j6  =  1/10.  The 
solid  line  again  is  the  exact  solution  and  the  dashed  line  is  the  3-region  result.  The 
peak  growth  rate  of  the  exact  solution  is  roughly  five  times  greater  than  that  of  the 
3-region  approximation  and  the  marginal  wavenumber  is  significantly  greater  than 
the  approximate  result.  We  emphasize  that  the  2-region  approximation  is  already 
completely  invalid  in  this  region,  giving  growth  rates  one  to  two  orders  of  magnitude 
lower  than  that  indicated  by  the  dashed-line  3-region  result.  Clearly,  the  3-region 
approximation  which  included  the  ion-intermediate  region  provided  significant  im¬ 
provement  over  the  conventional  approximations,  but  evidently  still  fell  far  short  of 
the  exact  solution.  The  shortcoming  of  the  3-region  approximation  can  be  traced  to 
the  fact  that  the  different  types  of  orbits,  crossing  and  non-crossing,  are  incorrectly 
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accounted  for  in  the  inner  and  outer  regions  and  that  the  eigenfunction  ij)  is  assumed 
to  be  constant  in  the  electron-inner  region  and  ion-intermediate  region.  These  are 
two  aspects  of  the  problem  that  cannot  be  easily  dealt  with  approximately  so  that 
an  integro-differential  approach  is  necessary.  In  Sec.  V,  we  will  discuss  these  issues 
in  greater  detail  in  formulating  a  simple  method  of  solving  the  integro-differential 
equation  while  correctly  taking  into  account  the  essential  physics. 

In  Fig.  •'5,  we  have  plotted  (solid  lines)  the  peak  growth  rate  versus  T^i/Tin 
for  several  values  of  Pi/S.  The  short  and  long  dashed  lines  represent  a  number 
of  approximations  and  will  be  discussed  later  (Sec.  V).  The  solid  curve  “a”  is  the 
case  Pi/ 8  =  1/2,  curve  “b”  corresponds  to  piJ8  =  1/10  and  curve  “c”  is  pi/ 8  = 
1/100.  Because  of  the  approximate  scaling  ')/ij^ci  ~  in  the  isotropic  limit^ 

the  growth  rate  is  greater  for  larger  pi/8.  However,  at  a  moderate  temperature 
anisotropy,  r,i/T,||  —  1  >  p,/^,  the  growth  rate  scales  as  an  inverse  power  of  pi/8. 
(The  approximate  theory  has  the  growth  rate  scaling  of  pi/ 8  in  this  regime.)  For 
small  values  of  pi/8  the  transition  between  these  two  regimes  of  anisotropy  occurs  at 
smaller  anisotropies,  so  that  the  growth  is  enhanced  more  quickly  as  the  anisotropy 
is  increased  in  this  regime  of  p,.  Thus,  curve  “a”  shows  the  least  variation  and  curve 
“c”  shows  the  most  variation  as  Tjx/Tlil  is  increased.  For  curve  “c”  the  growth  rate 
increases  by  nearly  four  orders  of  magnitude  as  Tix/T’,:||  increases  from  exact  isotropy 
to  l.rj. 

Figure  6  shows  the  peak  wavenumber  versus  r,x/r,||  for  the  examples  shown 
in  Fig.  5.  Again,  the  solid  curves  represent  the  exact  solution  and  the  short  and 
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long  dashed  lines  are  approximate  results  to  be  discussed  later.  For  the  present 
discussion,  the  important  point  is  that  the  wavenumber  corresponding  to  maximum 
growth  increases  with  increasing  favoring  shorter  wavelength  modes.  This 

effect  is  also  stronger  for  smaller  pJS.  The  results  of  the  3-region  approximation 
(short  dashed  lines)  are  in  qualitative  agreement  with  the  exact  result  but,  as  with 
the  dispersion  relation,  essential  quantitative  disagreement  develops  except  for  rel¬ 
atively  small  degrees  of  anisotropy. 

We  have  seen  that  the  collisionless  tearing  mode  in  a  neutral  sheet  requires 
solution  of  an  integro-differential  equation  except  for  limited  parameter  values 
(.5  ^  pijb  ^  .05  and  -  1  is  .5/o,/6).  However,  the  integro-differential 

treatment  is  highly  non-trivial  even  for  the  relatively  simple  system  given  here.  For 
more  complex  systems,  such  as  the  sheared  magnetic  field  geometry,  solution  of  the 
full  integro-differential  equation  may  be  prohibitive  in  terms  of  the  necessary  com¬ 
putation.  It  is,  therefore,  desirable  to  devise  a  method  which  provides  a  solution 
of  the  integro-differential  equation  without  sacrificing  accuracy  and  validity.  In  the 
following  sections,  we  will  describe  a  method  using  a  new  piecewise-straight-liiie 
approximation  of  crossing  orbits  in  evaluating  the  matrix  elements.  It  is  important 
to  point  out  that  this  method  still  solves  the  integro-differential  equation,  but  the 
time  required  for  computation  of  the  matrix  elements  can  be  reduced  significantly, 
allowing  for  solution  of  the  proljlem  in  complex  systems. 

In  order  to  understand  the  physical  basis  for  the  method,  we  will  use  the  3-region 
consideration  to  illustrate  the  essential  physics  required.  Then  we  will  apply  the 
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new  simplified  treatment,  to  the  bi-Maxwellian  tearing  mode  problem  solved  in  the 
preceding  sections  to  demonstrate  its  validity  for  a  wide  range  of  parameter  values. 
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IV  Review  of  the  “3-Region”  Approximation 


We  briefly  review  the  “3-region”  approximation  of  Chen  and  Paimadesso’,  which 
has  been  used  to  discuss  the  ion  anisotropic  tearing  mode.  In  this  approximation, 
the  sheet  is  divided  into  three  regions;  1.)  the  electron-inner  region,  |r|  <  within 
which  both  the  electrons  and  ions  are  modeled  by  free  trajectories,  2.)  the  ion- 
intermediate  region,  d,  '  |r|  '  d,,  where  the  electrons  are  magnetized  but  ions  still 
have  free  trajectories,  and  3.)  the  outer  region,  |c|  '  d,  where  both  species  are 

magnetized.  The  parameters  d,  and  de  are  chosen  to  be 

d,  =  .Jp,/(26). 

The  position  £  =  d^  is  the  turning  point  of  the  axis  crossing  electron  orbit  with 
p  =  0,  in  the  limit  <'<[  1,  and  d,  is  one  half  the  ttirning  point  of  the  orbit  with 

p  =  0. 

In  the  electron  inner  region  the  current  is  taken  to  be 

J\y  ~  Jad  +  • 

In  the  ion-intermediate  region,  we  have 

J\y  —  Jad  T  «/,  , 

and  in  the  outer  region 

Jly  —  J<id  +  +  J,  • 

Here,  the  adiabatic  current  is 

2T^5sech'(£),  (21) 
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For  each  species,  the  inner  region  current  is 


= 


and  the  outer  region  current  is 


J  =- 


47r^^  Tit  4-  Tie 


«2(0-(l  +  p^)(l+f2(0) 


/  1  \  sech"*!  i ) 

Had  we  included  tlie  density  profile  sech^(r)  in  writing  Eq.  (22),  the  resulting  matrix 
elements,  A’„„,  would  he  proportional  to  f^jdz  (?!>,„( £)sech^(i)/fj dr'  d>„(r').  Such  a 
non-symmetric  matrix  is  undesirable  on  both  physical  and  computational  grounds 
(the  eigenvalues  of  the  integro-differential  operator,  and  also  of  the  matrix  Gmr, 
must  be  reaF'.  The  lack  of  symmetry  arises  because  the  straight-line  model  orbits 
in  the  inner  region  are  incfuisistent  with  an  equilibrium  sech^(r)  pressure  profile. 
We  must  therefore  take  a  uniform  pressure  profile. 

However,  if  a  differential  equation  approximation  is  desired  (as  in  ref.  [l]),  we 
may  take  the  constant-i/'  approximation  for  the  integral  di'V’(i').  Then  there  is 
no  problem  with  retaining  the  physical  pressure  profile,  and  the  inner  region  current 


may  be  written  as 


^  r;)(' (24, 

Therefore,  if  we  model  the  inner  region  orbits  as  straight  lines,  we  must  either  take 
the  constant-V’  approximation,  or  we  must  use  a  constant  pressure  profile.  This 
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restriction  applies  because  the  inner  region  is  entirely  composed  of  crossing  orbits 
in  this  approximation.  When  we  discuss  the  piecewise-straight-line  approximation, 
which  does  not  have  an  inner  region  or  outer  region  as  such,  we  need  not  be  con¬ 
cerned  about  the  consistency  between  the  orbits  and  the  density  profile. 
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V  Features  of  the  Exact  Solution  and  Compari¬ 
son  to  the  3-Region  Approximation 


An  important  diagnostic  of  the  integro-differential  equation  is  the  effective  potential, 
defined  by 

V,fs{z)4iz)^ - 

c 

so  that  the  integro-differential  equation  takes  the  form  of  the  Schrodinger  equation 

d'^4' 


where  E  =  —k"^.  The  approximate  3-region  effective  potential  we  have  used  is  given 


■  Vo  -f  V- 

-f  v^< 

0  <  |£|  <  de 

K//(-)  —  ' 

v;,  +  V7 

+  V7 

de  <  \z\  <  di 

(25) 

,  Vo  +  v;> 

+  v;> 

d,'  1  1 1 

where 


V'  = 


Vo  -  -2  sech^(£), 


V.’ 


+  l)  sech^(£), 


Ti,  +  Ti 


V'  = 


2Tu 


n,  +  Tu 


tanh^(£)’ 

1  P 


and 


Vz’  = 


2Tu 


Ti,  -f  Ti, 


(26) 

(27) 

(28) 
(29) 


24 


This  potential  is  found  by  dividing  the  current  [Eqs.  (21), (24),  and  (23)]  by  V’(^)- 
In  Fig.  7,  we  show  the  effective  potential  for  a  solution  of  the  integro-differential 
equation  (solid  line)  together  with  the  a  the  3-region  approximation  to  the  effective 
potential  (dashed  line)  for  an  isotropic  example.  Although  the  overall  shapes  of  the 
exact  and  approximate  potentials  are  quite  similar,  important  differences  exist  both 
within  the  electron-inner  region  and  in  the  ion-intermediate  region.  For  isotropic 
plasmas,  we  will  find,  the  most  important  differences  between  the  exact  and  ap¬ 
proximate  potentials  are  due  to  the  response  of  electron  crossing  orbits  and  occur 
within  the  electron-inner  region.  The  differences  between  the  two  potentials  which 
exist  in  the  ion-intermediate  region  w’ill  be  found  to  be  of  less  importance  to  the 
isotropic  tearing  instability,  but  of  major  importance  to  the  anisotropic  instability. 

An  enlargement  of  Fig.  7,  which  shows  the  details  of  the  ion-intermediate  region 
potential,  is  shown  in  Fig.  8.  The  straight  long-dashed  line  is  the  level  E  —  —k^. 
In  the  exact  effective  potential,  there  are  additional  features  not  given  by  Eq.  (25). 
Specifically,  there  is  a  negative  spike  in  the  exact  potential,  located  just  outside  of 
the  electron  inner  region,  which  is  due  to  the  contribution  of  figure-8  electron  orbits. 
(Orbits  such  that  0  >  p  >  -y/h  make  a  figure  “8”  in  the  reference  frame  moving 
with  the  drift  speed  of  the  orbit.)  Such  orbits  are  not  fully  accounted  for  by  the 
straight-line  approximation.  Also,  the  ion  crossing  orbit  contribution  is  extended  to 
greater  values  of  £  than  in  the  approximate  potential.  (This  contribution  appears 
as  the  difference  between  the  solid  line  and  the  short  dashed  lane  from  £  =  d,  to 
£  ^  0.6.)  This  is  due  to  the  crossing  ions  which  have  energies  greater  than  the 
thermal  energy. 
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The  effective  potential  curve  for  an  approximation  may  be  used  to  assess  the 
inaccuracy  in  the  eigenfunction  solution.  The  area  between  the  approximated  and 
exact  potential  curves  relative  to  the  total  area  between  the  exact  potential  curve 
and  the  energy,  E  =  is  a  measure  of  the  inaccuracy  incurred  in  using  the 

particular  approximation.  In  addition,  the  inaccuracy  attributable  to  a  particular 
type  of  orbits  can  be  evaluated  by  calculating  the  area  difference  due  to  the  or¬ 
bits  in  question.  Features  in  Fig.  8,  specifically  the  negative  spike  due  to  electron 
figure-8  orbits  and  the  contribution  from  the  high  energy  ion  orbits  are  found  to  be 
unimportant  in  the  isotropic  case. 

Of  much  greater  importance  is  the  difference  between  the  exact  potential  and 
the  approximate  potential  within  the  electron  inner  region,  as  was  shown  in  Fig.  7. 
The  constant- v'’  approximation  is  probably  the  origin  of  the  large  difference  between 
the  potential  of  the  3-region  approximation  and  the  exact  potential  in  the  electron- 
inner  region.  In  reality,  since  the  potential  is  much  larger  than  than  the  energy 
the  eigenfunction  is  strongly  evanescent  in  the  electron-inner  region  and  not  at 
all  constant.  The  physical  interpretation  is  this:  Rather  than  a  uniform  perturbed 
current  density  forming  throughout  the  electron-inner  region,  actually  strong  “skin” 
currents  form  at  the  edge  of  the  inner  region  as  a  response  to  the  induced  electric 
field  Ey  =  -yi'/c.  However,  a  certain  amount  of  total  current  is  required  within  the 
electron  inner  region  to  maintain  the  magnetic  field  perturbation,  and  this  total 
current  is  now  seen  to  be  carried  by  the  thin  “skin”  current,  rather  than  by  a 
uniform  current  density.  Thus  the  induced  electric  field  must  be  larger  than  the 
constant-V’  approximation  would  lead  us  to  believe,  implying  a  larger  growth  rate. 
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Therefore,  we  see  that  the  constant-^’  approximation  under  estimates  the  growth 
rate  because  it  does  not  account  for  the  current  density  gradient  in  the  inner  region. 

The  growth  rate  for  the  isotropic  case,  pi/8  =  1/10,  as  a  function  of  k  is  given 
in  Fig.  3.  Although  the  3-region  approximation  (dashed  line)  under  estimates  the 
growth  rate,  as  has  been  predicted,  the  approximation  does  much  better  here  than 
in  some  cases  we  will  discuss  presently. 

In  the  anisotropic  case,  shown  in  Fig.  4,  the  agreement  of  the  3-region  approx¬ 
imation  with  the  exact  solution  is  not  nearly  so  good.  Here,  the  peak  growth  of 
the  exact  solution  is  about  four  times  the  peak  of  the  3-region  approximation  and 
unlike  in  the  isotropic  case,  the  3-region  approximation  now  over  estimates  the 
wavenumber  corresponding  to  the  peak  of  the  growth  rate. 

The  effective  potential  for  the  fastest  growth  eigenfunction  in  Fig.  4  is  shown  in 
Fig.  9.  The  minimum  value  of  the  effective  potential  in  the  ion-intermediate  region 
is  roughly  given  by  the  3-region  approximation,  which  for  these  parameters  scales 
with  1  —  T_l/T',||,  however,  important  differences  exist  in  the  shapes  of  these  two 
potentials.  For  p,/6  small,  as  in  this  case,  the  contribution  from  non-crossing  orbits 
is  small.  The  contribution  to  the  potential  from  the  crossing  orbits  is  therefore  the 
difference  between  and  the  adiabatic  potential  Fq  =  sech^(£);  in  the  3-region 
approximation,  this  difference  is  given  by  .  It  is  important  to  note  that  in  the 
anisotropic  case  the  contribution  from  the  ion  crossing  orbits  has  changed  sign, 
and  again  extends  to  larger  values  of  z  than  the  3-region  approximation  indicates. 
Further,  the  contribution  has  increa.sed  in  absolute  value,  so  that  the  area  between 
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the  exact  and  approximate  potential  in  the  ion-intermediate  region  is  comparable 
to  the  total  area  between  the  exact  potential  and  the  energy.  Thus,  the  differences 
between  the  exact  and  approximate  potentials  due  to  the  ion  crossing  orbits  domi¬ 
nate,  rather  than  differences  due  to  electron  crossing  orbits  which  dominate  in  the 
isotropic  case. 

The  effect  of  increasing  the  ion  temperature  anisotropy  on  the  0  eigenfunction 
is  shown  in  Fig.  10  (these  are  the  eigenfunctions  corresponding  to  the  maximum 
growth  rate).  As  the  anisotropy  is  increased  the  eigenmode  appears  to  becomes 
more  strongly  trapped  in  the  ion-intermediate  region  potential  well,  which  becomes 
deeper.  Because  the  eigenfunction  becomes  more  trapped  in  the  ion-intermediate 
region,  constant-0  approximation  breaks  down  inside  of  |£|  =  d,. 

To  summarize  the  comparison  between  the  isotropic  and  anisotropic  case,  in 
the  anisotropic  case  the  ion  contribution  to  the  effective  potential  (K//  —  ho)  in 
the  ion-intermediate  region  has  changed  sign  and  the  difference  between  the  ion- 
intermediate  region  potential  well  and  the  adiabatic  potential  well  is  greater  than 
in  the  isotropic  case.  Thus  the  difference  between  the  potential  of  the  three  region 
approximation  and  the  potential  of  the  the  exact  solution  in  the  ion-intermediate 
region  will  become  more  important  as  the  anisotropy  is  increased.  This  difference 
is  primarily  due  to  the  inclusion  into  the  exact  potential  of  the  contribution  of 
high  energy  crossing  orbits.  Further,  the  constant-0  approximation  of  the  3-i'egion 
approximation  will  be  less  valid  as  the  anisotropy  is  increased. 

Returning  to  Figs.  5  and  6,  we  can  see  that  as  the  temperature  anisotropy  is 
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increased,  the  3-region  approximation  (dashed  line)  disagrees  more  with  the  exact 
solution  (solid  line).  Also,  as  Pi/S  is  decreased,  the  disagreement  of  the  3-region 
approximation  becomes  greater.  As  we  can  see  from  the  effective  potential,  shown 
in  Fig.  11,  as  pi/ 8  is  decreased,  the  contribution  to  the  effective  potential  from 
crossing  orbits,  especially  the  electron  crossing  orbits,  is  increased.  Thus,  as  pi/b 
is  decreased,  the  role  of  the  electron  inner  region  “skin”  currents  discussed  earlier 
becomes  greater  and  the  constant-V’  approximation  becomes  more  invalid.  Thus  we 
have  the  rather  surprising  result  that  as  /?,/6  is  decreased,  the  constant-^’  approxi¬ 
mation  is  less  valid. 

The  dashed  line,  which  agrees  with  the  solid  line  much  better,  especially  in  the 
large  p,  (pt/^  =  1/2)  and  in  the  small  p,  (p,/6  =  1/100)  cases,  gives  the  frequency 
solution  of  our  piecewise-straight-line  approximation. 

A  A  piecewise-straight-line  approximation 

We  have  solved  the  integro-differential  equation  in  the  neutral  sheet  geometry  with¬ 
out  approximating  the  particle  orbits  or  the  perturbed  potential  eigenfunction,  and 
the  results  show  that  in  much  of  parameter  space  either  the  constant-^’  approxi¬ 
mation  is  invalid  or  the  approximate  orbits  as  in  the  3-region  approximation  are 
not  valid.  Computation  of  the  A’m„  matrix  elements,  however,  is  time  consuming. 
In  fact,  if  the  matrix  elements  depended  upon  more  than  a  single  parameter,  as 
they  fortunately  do  in  this  geometry,  the  problem  would  not  be  practical  to  solve 
in  this  way.  In  more  complicated  geometries,  such  as  the  sheared  field  geometry, 
the  matrix  elements  would  not  only  be  more  complicated  to  compute,  but  they 
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would  have  to  be  recomputed  more  frequently  to  perform  a  parameter  search.  In 
the  sheared  field  geometry,  for  example,  the  matrix  elements  would  also  depend 
upon  the  wavenumber  as  well  as  the  dependence  upon  the  asymptotic  gyroradius. 

The  essentials  of  the  full  integro-differential  equation  is  that  no  form  for  the  per¬ 
turbed  potential  eigenfunction  is  assumed,  such  as  the  constant-V’  approximation, 
and  that  orbits  are  treated  as  crossing  or  non-crossing  depending  upon  the  values  of 
the  azimuthal  momentum  and  the  energy,  rather  than  the  position  i.  The  adverse 
impact  of  approximating  crossing  orbits  by  straight  lines  seems  to  be  smaller. 

Therefore,  we  have  approximated  the  orbit  integration  over  a  basis  function  by 
assuming  that  the  velocity  is  roughly  constant  over  that  basis  function 


dz'  4>^[z') 


provided  the  orbit  is  a  crossing  orbit.  Here,  and  z^^x  are  the  minimum  and 

maximum  extents  of  the  chapeau  basis  function  centered  at  If  the  turning 
point  occurs  within  the  basis  function,  however,  the  contribution  from  that  orbit 
is  discarded.  It  would  also  be  straightforward  to  include  the  contribution  from  the 
turning  point  in  some  approximate  way.  The  exact  contribution  from  the  electrons, 
and  the  exact  contribution  from  the  non-crossing  ion  orbits  are  retained  in  this 
approximation. 

As  we  may  see  from  Figs.  5  and  6,  the  piecewise-straight-line  approximation 
agrees  very  well  in  both  the  small  and  the  large  pi/h  cases  with  the  exact  solution. 
Even  in  the  intermediate,  pi/b  =  1/10  case,  it  is  much  better  than  the  .3-region 
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approximation.  The  inclusion  of  the  turning  point  contributions  would  make  the 
agreement  of  the  piecewise-straight-line  approximation  much  better  in  all  cases. 
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VI  conclusions 


We  have  calculated  the  exact  growth  rate,  7(fc),  eigenfunctions,  and  effective  poten¬ 
tials  for  the  collisionless  tearing  mode  in  an  anisotropic  bi-Maxwellian  neutral  sheet 
by  solving  the  integro-differential  equation  and  using  the  exact  particle  orbits.  The 
results  are  not  qualitatively  different  than  those  reported  earlier  in  the  literature', 
although  the  best  previous  approximation,  the  3-region  approximation,  under  esti¬ 
mates  the  growth  rate  by  a  factor  of  two  to  five.  We  find  that  the  presence  of  a 
favorable  anisotropy  (T.x/Tin)  greatly  increases  the  growth  rate  and  decreases  the 
characteristic  wavelength,  especially  in  the  case  of  pi/S  <C  1. 

When  comparing  the  numerical  integro-differential  equation  solution  to  the  3- 
region  approximation  for  the  isotropic  tearing  mode  (see  figs.  12  and  13)  we  found 
that  the  3-region  approximation  breaks  down  in  the  small  p,  regime.  In  this  regime, 
the  increase  in  the  difference  between  the  height  of  the  effective  potential  in  the 
electron  inner  region,  and  the  adiabatic  potential  means  that  the  eigenfunction  is 
strongly  evanescent  inside  of  the  electron  inner  region,  and  constant-V’  approxima¬ 
tion  of  the  3-region  approximation  breaks  down. 

When  an  ion  anisotropy  was  included,  the  difference  between  the  three  region 
approximation  and  the  numerical  solution  of  the  integro-differential  equation  also 
increased.  This  is  due  primarily  to  the  extension  of  the  contribution  of  crossing 
ions  to  the  potential  for  positions  z  outside  of  the  ion-intermediate  region  which  is 
due  to  high  energy  crossing  ions,and  to  the  break  down  of  the  local  approximation 
in  the  ion  intermediate  region.  In  particular,  for  small  pi,  the  eigenmode  becomes 
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deeply  trapped  in  the  ion-intermediate  region  as  the  ion  anisotropy  is  increased,  and 
therefore  an  approximation  which  assumes  V’  is  constant  within  the  ion-intermediate 
region,  such  as  the  3-region  approximation,  is  contradictory. 

We  have  also  used  a  piecewise-straight-line  approximation  to  solve  the  integro- 
differential  equation  in  a  simplified  manner.  This  approximation  allows  quick  eval¬ 
uation  of  the  matrix  elements,  facilitating  the  solution  of  the  integro-differential 
equation,  although  it  retains  the  essential  features  of  the  integro-differential  treat¬ 
ment.  In  particular,  this  piecewise-straight-line  approximation  correctly  classifies 
the  orbits  as  crossing  or  non-crossing  based  on  their  value  of  momentum  and  energy 
rather  than  upon  position,  and  assumes  no  particular  form  for  the  eigenfunction. 
This  should  be  contrasted  with  the  conventional  straight-line  and  constant-V^  ap¬ 
proximation  in  which  all  orbits  within  a  given  distance  from  the  null  plane  are 
treated  as  straight-line  orbits  without  distinguishing  between  crossing  and  non¬ 
crossing  orbits  and  in  which  i/’  is  assumed  to  be  constant  near  the  null  plane. 

For  our  piecewise-straight-line  approximation,  we  assumed  that  the  velocity  is 
constant  over  one  basis  function  for  ion  crossing  orbits.  When  a  turning  point 
occurred  within  a  basis  function,  however,  we  ignored  the  contribution  from  that 
orbit.  A  further  refinement  to  this  approximation  would  be  to  include  the  turning 
point  in  some  approximate  fashion.  In  the  large  p,  limit,  our  approximation  becomes 
exact  since  these  crossing  orbits  are  much  larger  than  a  basis  function.  Because  in 
the  small  p,  limit,  the  region  that  the  crossing  orbits  occupies  is  small,  the  details 
of  particle  orbits  are  not  so  important.  Thus,  in  the  small  pi  limit  the  key  aspect  of 
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the  piecewise-straight-line  approximation  is  that  is  does  not  assume  any  particular 
form  for  the  t/’  eigenfunction. 
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Figure  1.  The  region  of  integration  for  some  non-crossing  orbits.  The  gray  shaded  region 
is  the  region  of  integration  for  the  non-crossing  orbit  with  p  =  —y/h,  while 
the  cross-iiatched  area  is  for  an  orbit  wi  th  p  <  —y/h. 
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ure  2.  The  region  of  integration  for  some  crossing  orbits.  The  gray  region  is  the  area 
of  integration  for  the  crossing  orbit  with  the  largest  turning  point,  p  =  -\/h. 
The  cross-hatched  region  is  for  an  orbit  with  smaller  turning  points,  i.e.  p  > 

-Vh. 
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Figure  5.  The  peak  growth  rate  as  a  function  of  ion  anisotropy,  “a”,  pi/b  =  1/2.  “b”, 
Pi/S  —  1/10.  “c”,  pi/f)  —  1/100.  The  short  dashes  are  the  solution  with 
the  3-region  approximation,  the  long  dashes  are  found  using  out  piecewise- 
straight-line  approximation,  described  in  Sec.  A. 


Figure  8.  An  enlargement  of  Fig.  3,  showing  the  ion-intermediate  region  potential  well. 
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Figure  10.  Eigenfunctions  with  pijb  —  1/10.  “a”;  7’,x/T',||  =  1,  fc  =  0.31  and 

1..58  X  10-^  “b”;  r,x/r,||  -  1.0.5,  k  =  0.52  and  =  3.74  x  lO"'*.  “ 

r,x/T,||  =  1.2,  k  =  0.86  and  7/u;e,  =  1.96  x  10"^  “d”;  r,x/r,||  =  1.5,  k  =  1, 
and7/a;,.  =  1.44x10-^  “b”;  r,;x/7;||  =  1.05,  k  =  0.52  and  7/u>e,  =  3.74x10 


Figure  11.  The  effective  potential  well  with  pijb  =  1/100. Tix/Tin  =  1,  ^  =  ..5,  and 
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